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The energy of a displaced magnetic vortex in a cylindrical 
particle made of isotropic ferromagnetic material (magnetic 
dot) is calculated taking into account the magnetic dipolar 
and the exchange interactions. Under the simplifying assump- 
tion of small dot thickness the closed-form expressions for the 
dot energy is written in a non-perturbative way as a function 
of the coordinate of the vortex center. Then, the process of 
losing the stability of the vortex under the influence of the ex- 
ternally applied magnetic field is considered. The field desta- 
bilizing the vortex as well as the field when the vortex energy 
is equal to the energy of a uniformly magnetized state are 
calculated and presented as a function of dot geometry. The 
results (containing no adjustable parameters) are compared 
to the recent experiment and are in good agreement. 
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Introduction. The unusual magnetic properties of 
sub-micron cylindrical magnetic particles and their peri- 
odical two-dimensional arrays drawn much attention be- 
cause of their possible potential as a magnetic storage as 
well as an interesting jmadel system for studying of the 
magnetization reversaJJTa. The applications of such pat- 
terned magnetic films for magnetic information storage 
(MRAM, for instance) are promising. 

In this work the isolated magnetic particles made 
of isotropic (soft) magnetic material shaped as circular 
cylinders (referred hereafter as "dots") are considered. 
The magnetic properties of the dots are governed by the 
magnetic dipolar and the exchange interactions. When 
the external magnetic field is absent, three parameters 
of the dimension of length define completely the mag- 
netic structure of a dot. These parameters are: the dot 
radius R, the dot thickness L and the exchange length 
Le = y C/Mg, where C is the exchange constant and 
M$ is the saturation magnetization of the material. The 
phase diagram of small L/Le < 4, R/Le < 4 dots in 
no applied magnetic field showing the range of dot pa- 
rameters where the magnetic vortex (curl) is rthje ground 
state was calculated by Usov and co-workersOu. Vortex 
states—wese indeed observed in polycrystalline Ccu and 
FeA^fl-y cylindrical dots. 

So far there is no theory of magnetization reversal 
of a dot in a vortex state when the magnetic field 
h = H/4ttMs applied in the direction perpendicular to 
the cylinder axis. This process taking place via nucle- 
ation, displacement and annihilation of a single magnetic 
vortextil is considered in this work. There are three char- 
acteristic fields describing this process: the "vortex" nu- 



cleation h n and annihilation /i an critical fields (for tran- 
sition to/from the curling state from/to the saturated 
state) and also the field h eq , corresponding to the vor- 
tex state having the same energy as the uniform one. 
The size dependent fields of vortex equilibrium h cq and 
annihilation h an are calculated in this paper. The field 
of vortex nucleation h n is determined by the imperfec- 
tions of the dot shape (and/or material) and was exper- 
imentallvpipund to have nearly no correlation to the dot 
geometrjO. The field h cq enters directly the phase di- 
agram in coordinates h, L/Le, R/Le whose zero field 
partfl was already presented. 

A realistic description of the non-uniform magnetiza- 
tion distribution in an isolated dot at zero applied field 
was already considered: first, with the help of guessed 
trial functionsEj and then using the variational principled. 
These approaches heavily rely on the cylindrical coordi- 
nate system and are hard to generalize for the case when 
the vortex is deformed by the in-plane magnetic field. 
Recently, another variational method to construct the 
magnetization distributions in small exchange-dominated 
cylinders (not limited to circular ones) was presentedO. 
It allows to obtain the structure of a displaced vortex, 
which is a starting point for the calculations done in this 
work. 

Equilibrium parameters and stability of displaced vor- 
tex. In a thin ferromagnetic cylinder, with the thickness 
less or of the order of the exchange length, the magneti- 
zation distribution is constant along the thickness and is 
essentially two-dimensional. In the Cartesian coordinate 
system with coordinates {A, Y, Z}, chosen in such a way 
that the axis 0Z is parallel to the cylinder axis, the nor- 
malized magnetization vector field m(r) = M(f)/Ms, 
\m\ = 1 is independent on the coordinate Z, so that 
r = {A, Y}. In this case the magnetization distribution 
corresponding to the minimum of the exchange energy 
functional 
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where the components of the magnetization vector are 
expressed as m x + rm y — 2w/(l + ww) and m z — 
(1 — ww) /(l + ww) through the complex function w(z, z) 
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of the complex variable z — X/R + iY/R, the horizon- 
tal line over a variable denotes the complex conjugation 
z = Xj R — iY/R, the straight brackets stand for the ab- 
solute value of a complex number \z\ 2 = z~z. The expres- 
sion (|^) gives the magnetization distribution minimizing 
the functional ([[]) for all values of a scalar parameter c, 
related to the radius of the vortex core and a complex 
parameter a = Xc/R + lYc/R, describing the position 
of the vortex center. The_case of a = corresponds to 
the Usov's vortex solution^ which satisfies the conditions 
(fh-n)s — and (dm/dn)s = on the side dot boundary. 
The conformal mapping z (z — a)/(l — az) usedlla to 
obtain (2) holds the boundary condition (dm/dn)s = 
unchanged, whereas the side magnetic charges (m • n)s 
appear. If the magnetic field is applied along Y axis 
it is clear from the symmetry that the vortex will be 
displaced parallel to the X axis, which allows to as- 
sume that the parameter a is real. The expression (Q) 
yields correct limiting cases a = (centered vortex) and 
a = 1 (single-domain or saturated dot) and is applicable 
at small enough R and L. The example distribution at 
an intermediate value of a is given in Fig. [L| To find par- 
ticular values of the parameters a and c corresponding to 
the given dot geometry, material and applied magnetic 
field it is necessary to include into consideration the mag- 
netostatic and Zeeman energies. 

The magnetostatic energy has two surface (on the dot 
faces and sides) and a negligible (for small dot thickness) 
volume contributions. To simplify the further consid- 
eration we shall assume that the dot thickness is small 
enough to neglect the volume magnetostatic contribu- 
tion. 

In the further calculations two polar coordinate sys- 
tems with the same notation for coordinates r, <J> will 
be used. One of the coordinate systems will be centered 
in the dot center and have the dimensionless polar ra- 
dius p = r/R. The other coordinate system centered 
at the vortex core center is also polar, because the vor- 
tex core (the region of the dot with non-zero m z ) has 
exactly the circular boundary for all values of a. The 
center of the vortex core (note the difference with the 
vortex center where m z = 1) is situated at the coordi- 



nates Xy = a(l — & 
its radius Ry = c(l 
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= for real a, and 
The dimensionless 
polar radius for the coordinate system centered at the 
vortex core is p = r/Ry. 

The density of magnetic charges produced by the dis- 
tribution (0) on the faces of the dot in the polar coordi- 
nate system centered at the vortex core is 



m z =cr(p,t 
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with 5 = 2 arctan(ac) . The magnetostatic energy of these 
charges normalized to 4%MgV (V = LttR 2 is the dot 
volume, all energies further in the text will be given in 
this normalization) is 
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where 
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the integration runs from to 1 in p\, p 2 and from to 
2tt in cf>i, cf>2- Assuming that the parameter ac is small 
(which is justified by our numerical calculations, also see 
the caption to Fig. [l]) the last expression in zero-th order 
in ac is 
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where Gu(x) = G(x,S — > 0) is the magnetostatic func- 
tion of non-displaced vortex corea and g = L/R. While 
this expression is approximate it captures the details of 
the exact one (|]) well, and has error not more than 1% 
as we have checked numerically. 

The surface magnetic charges on the sides of the dot 



(mri)s = c(0) = 2asin(/> 
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which is written in the cylindrical coordinate system r, (f> 
centered at the dot center. The energy of these charges 
can be expressed as 
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where the integration runs from to g in z\ ,Z2 and from 
to 2-7T in 0i, <p2- The last expression can be simplified 
using the summation theorem for the Bessel's functions 
into 
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where J p (x) is the Bessel's function of the first kind and 
order p,. Minimizing the total energy we have kept only 
one term in the above series, which was sufficient to have 
the error less than 3%. 

The total magnetization of the dot (having only one 
non-zero component along Y axis) , which enters the Zee- 
man energy is 



((1 + a 2 )pcoscj)— a(l + p 2 ))pdpd(j) 



\J 'a 2 + p 2 — 2ap cos (f>\/l + a 2 p 2 — 2ap cos <fi 



(11) 



where the integration from to 1 in p and from to 27r 
in 4> is performed in the coordinate system centered on 
the dot. This integral can be approximated by (m y ) = 
a(l5a 4 + Ma 2 - 193)/144 with the error less than 1% 
according to our numerical calculations. 

The ex&hange energy can be evaluated using its 
expression!! 3 ] directly through the function w(z, ~z), or, al- 



tiy 

ternatively, using (|]) as 



(12) 
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which is independent on a. 

The total energy density of the particle is e = e ox + 
e MS + e Ms ~ h( m y) an d should be minimized with re- 
spect to the parameters a and c. These parameters are 
coupled through the e^fj term which depends on both 
of them. However, if we use q = c(l — a 2 )/ g instead of c 
the parameters a and q become decoupled and the energy 
can be minimized with respect to each independently. 

The only terms depending on q are eEx + ^ms = 
X 2 g 2 (2 - \ogqg) + q 3 g 2 Q(l/q) + const, where Q{x) = 
Gu(0) — Gjj(x), A = Le/{L\^tt). The equilibrium value 
of q, then, satisfies the following equation 



q fq ^(Vff)] = A ' 



(13) 



which is exactly the equatioroQ for the equilibrium radius 
of the vortex core at zero field. 

The rest of terms in the total energy ei(a, h) — effg — 
h(m y ) + X 2 g 2 log(l — a 2 ) depend only on a, the last term 
comes from the exchange energy after changing the inde- 
pendent variable from c to q. 

The vortex loses its stability when at a certain value 
of the external field h = h aa the energy minimum it is in 
turns to a saddle point or a maximum. This leads to the 
following system of equations for et an , h an : 
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These equations were solved numerically for Lr&JL = 
18 nm/15 Tim (corresponding to the experimentO) and 
the resulting dependence of the vortex annihilation field 
h an on the reduced dot radius 1/g = R/L is shown in 
Fig. ||. The resulting-/7 an (l/(?) curve agrees well with the 
experimental pointsO. 

The field h cq when the energy of the displaced vortex 
in the external magnetic field is equal to the energy of 
the uniformly magnetized state was determined from the 
equations: 



<9ei(a e q,/ieq)/<9a eq = 
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and also plotted in Fig. This expression takes into 
account the fact that the exchange and the face magne- 
tostatic energies are equal to zero for the single-domain 
(saturated in-plane) dot. 

Discussion. The calculations of the annihilation field 
explain reasonably the experiment over a wide range of 
1/g. The differences between the measured and calcu- 
lated h an curves at small R/L < 10 can be the result of 
increasing role of thermal fluctuations in small particles. 
The peak on the calculated dependencies of both h an and 
h cq on R (see Fig. |^) due to their decrease with the cylin- 
der cadius at small R is not present in the experimental 
datatU and needs an explanation. 

It can be explained as the result of inherent (not requir- 
ing the applied field) instability of vortex in the cylinders 
of small radius. There are two kinds of forces acting on 
the vortex in a cylinder, the stabilizing (due to the energy 
terms increasing as the vortex comes off the dot center) 



and destabilizing ones. In the case considered in this pa- 
per the main stabilizing force is produced by the energy 
of magnetic charges on the cylinder side. The rate of the 
increase of this energy term (||) with a (the actual force) 
decreases with the dot radius R at a fixed L. This de- 
crease is related to the decrease in the possible amount 
of magnetic charges forming (the area of side decreases 
with R) and also to the decreasing distance between the 
positive and negative ones, whose negative interaction 
energy strengthens as a result. 

The destabilizing force is due to the energy of the mag- 
netic charges on the cylinder faces (|4|), which is positive 
for the centered vortex and is zero when the vortex is 
completely out of the cylinder. This force is independent 
on the dot radius as long as the vortex core fits com- 
pletely within the particle. 

Therefore, as the radius of the cylinder decreases there 
is a point (called the stability radius Rs) where the desta- 
bilizing force starts to outweigh the stabilizing one and 
the vortex becomes unstable in such a cylinder. If the 
cylinder radius is smaller than Rs no magnetic field is 
required to destabilize the vortex and /i all = 0. This 
means that the field h an decreases at small cylinder radii 
to become zero at R = Rs- 

The similar decrease of h cq happens because the 
smaller is the cylinder the smaller is the energy of the 
uniformly magnetized state in it and it is closer to the 
energy of the cylinder with a vortex. Until, at the so- 
called single domain radius Rc these energies become 
equal with no applied field. This means at small radii 
h eq also decreases with the radius until it is exactly zero 
at R = R c . 

The values of both Rs and Rc have sense only when 
they are larger then the vortex core radius Ry. This 
is not the case for the cylinder thickness L/Le = 15/18 
shown in Fig. g, where Rs < Rv)- Our calculations show 
that at smaller dot thicknesses the case when Rs , Rc > 
Rv can be realized and h an , h oq curves start from zero 
at small radii. 

As we have mentioned in the introduction, the de- 
pendence n of Rc/Le on L/Le was already reported in 
literature^, its calculation does not require a model of 
displaced vortex. The coincidence of the cylinder radius 
at which h eq turns to zero with Rc was one of the checks 
we did to verify the results. 

The approximations we performed were focused on cal- 
culating the dot energy at arbitrarily large vortex dis- 
placements, that is, not to perform a perturbative expan- 
sion limiting us to a particular range of displacements. 
However, for the calculation of Rs corresponding to the 
vortex displacement mode given by Eq. (|J) the expansion 
at a = up to second order is enough. The approxima- 
tions done here, especially Eq. (^]), are, strictly speaking, 
not valid up to the second order in a at a = 0. Therefore, 
the stability radius resulting from this calculation should 
be taken with caution. The mentioned expansion is easy 
to perform and we plan to present the precise estimate 
of Rs for the mode (||) in one of our forthcoming papers. 

Some correlation of the vortex annihilation field with 
the dot average in-plane demagnetization factor in satu- 
rated state was already notecfl. That can be explained as 



3 



follows. The vortex shift a in the external magnetic field 
is determined by competition of the side magnetostatic 
and Zeeman energies. The main contribution in Eq. (||) 
to the magnetostatic energy from the dot boundary is 
proportional to I\(g) being exactly the single-domain dot 
in-plane average demagnetization factorial. But this is 
correct only for large enough dot radius (considerably 
above Rs, Rc) when the face magnetostatic energy and 
exchange energy can be neglected. 

Conclusions. In summary, the size-dependent vortex 
equilibrium and annihilation fields have been studied the- 
oretically for circular cylindrical magnetic dots. Simple 
analytical approximation suggested for the dot magne- 
tization distribution allowed us to consider the transi- 
tion from the vortex to single-domain dot magnetization 
state under applied magnetic field. The calculation re- 
sults agree well with the experimental data within the 
limit of weak inter-dot coupling. Further studies of static 
and dynamic properties of magnetization in isolated and 
magnetostatically coupled dots are in progress. 
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FIG. 1. The magnetization distribution ([|) with the equi- 
librium values a » 0.14 and c « 0.092 corresponding to 
a permalloy (Le = 18 nm) cylinder with R = 150 nm, 
L = 15 nm in the magnetic field h = H/{A-kMs) = —0.025 
applied parallel to the cylinder OF axis. The out-of-plane 
magnetization component on the line Y = is shown on the 
upper plot. 
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FIG. 2. Characteristic fields of the permalloy Le = 18 nm 
dots of the thickness L — 15 nm as a function of their padius. 
The experimental points are taken from the literature!! 1 !. 
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